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We consider corrections to the underlying cosmology due to secondary contributions from weak 
gravitational lensing, the integrated Sachs- Wolfe effect, and the Sunyaev-Zel'dovich effect contained 
in the trispectrum. We incorporate these additional contributions to the covariance of a binned 
angular power spectrum of temperature anisotropies in the analysis of current and prospective data 
sets. Although recent experiments such as ACBAR and CBI are not particularly sensitive to these 
additional non-Gaussian effects, the interpretation of Planck and CMBPol anisotropy spectra will 
require an accounting of non-Gaussian covariance leading to a degradation in cosmological parameter 
estimates by up to 20% and 30%, respectively. 
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I. INTRODUCTION 

The cosmic microwave background's (CMB) ability to 
constrain the cosmological parameters has driven a large 
variety of CMB experiments. The high quality mea- 
surements of the temperature and polarization of the 
CMB are compatible with a flat universe with nearly 
scale invariant fluctuations. The sensitivity to small an- 
gular scales of recent experiments such as WMAP5 
ACBAR 0] and CBI give us a better understanding 
of anisotropies from local large scale structure (LSS) and 
from non-linear effects such as weak gravitational lens- 
ing. The non-Gaussianity of these fluctuations can have 
an effect in the estimation of the cosmological param- 
eters, especially for future experiments. However, pos- 
sible non-Gaussian effects have not been considered in 
previous analyses, under the assumption that they are 
negligible. 

Non-Gaussianities show up as contributions to the 
four-point correlation function, or trispectrum in Fourier 
space, of the CMB temperature fluctuations H, [H]. The 
four point correlations quantify the sample variance and 
covariance of the two point correlation of power spectrum 
measurements 0, 0] ■ To adequately understand the sta- 
tistical measurements of CMB anisotropy fluctuations a 
proper understanding of the four point contributions is 
needed. Given the high precision level of cosmological pa- 
rameter measurements by current and future surveys, a 
careful consideration must be attached to understanding 
the presence of non-Gaussian signals at the four point 
level. The goal of this paper is to understand to what 
extent these non-Gaussianities affect constraints on the 
cosmological parameters 

As explored in previous work (e.g. (8l-[TT|V one of 
the most important non-linear contributions to CMB 
temperature fluctuations results from weak gravitational 
lensing. Contributions to the four point level arise from 
the non-linear mode-coupling nature of the lensing effect, 
as well as correlations between weak lensing angular de- 
flections and secondary effects that trace the same large 
scale structure. The trispectrum due to lensing alone is 



studied in [5|, and further considered under an all-sky 
formulation in Q • 

In our analysis we include contributions to the trispec- 
trum from three well known secondary sources. First, 
we consider the contribution of the trispectrum to the 
power spectrum covariance due to lensing. We also con- 
sider contributions from the integrated Sachs- Wolfe effect 
(ISW; [11) and the Sunyaev-Zel'dovich effect (SZ; [111), 
both which cross-correlate with lensing to produce non- 
Gaussianity. Each of these effects becomes non-negligible 
on small angular scales and must be taken into account 
for a more complete treatment of cosmological parameter 
constraints. 

In a situation without instrumental noise, the lensing 
contribution to the trispectrum can increase the values 
along the diagonal of the covariance matrix on the or- 
der of about 10% The SZ contributions are greater 



than the lensing by as much as an order of magnitude. 
However, with noise added follows the possibility that 
these effects on the trispectrum become negligible. Most 
experiments thus far have been too noisy at small angu- 
lar scales to observe these non-Gaussian effects on the 
cosmological parameters. We explore how firmly this 
holds for the ACBAR and CBI data sets when combined 
with WMAP five-year data as well for future Planck and 
CMBPol measurements. 

Moreover, we consider the weak lensing scaling param- 
eter Al in our analysis. This parameter scales the power 
spectrum of the CMB lensing potential Cf"^ such that 
Al = Q corresponds to an unlensed scenario, whereas 
Al — 1 renders the expected lensed result. It was first 
reported that in an analysis of the ACBAR data set this 
parameter is inc omp atible with the expected value of 
Al ~ 1 at 2.5(7 [Taj. ACBAR has since updated their 
data set with a new calibration, and the ACBAR team 
finds that the amplitude is consistent with the expected 
value of unity [l6|. However, given that this parameter 
is sensitive to physics at high multipoles, such as that 
encapsulated in the trispectrum, we also include A^ as 
an extra parameter in our calculations. 

Our analysis of the ACBAR and CBI data sets, com- 
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bined with WMAP 5-year data, show that secondary 
contributions in the trispectrum have a negUgible ef- 
fect on the uncertainties of the cosmological parameters. 
At most the constraints on the cosmological parameters 
change by no more than 10%. We find that these modifi- 
cations are negligible due to the increase in noise at large 
multipolcs (as shown in Fig. [S|). 

However, in a mock Planck data set with realistic noise 
contributions, the non-Gaussian effects coming from the 
trispectrum become more important, with error con- 
straint degradations up to 20% (as shown in Table ITV]) . 
For the case of CMBPol this rises to < 30%. For this 
reason, future experiments with at least Planck-level sen- 
sitivity at high I should consider non-Gaussianities from 
the trispectrum in order to have proper error bar esti- 
mates. 

Lastly, we show Al to be consistent with the expected 
value of unity, primarily due to the same reasons that 
Ref. [l^ found it to be consistent with the theoretically 
expected value. If we had used the original ACBAR 
dataset, prior to their most recent update, we find that 
Al is inconsistent with unity at the same confidence as 
found by Ref. [l^. We also find that this parameter 
remains consistent with unity when both trispectrum ef- 
fects and the running dus/dlnk of the spectral index Ug 
are included in the analysis. 

In Sjlll we review the calculation for the trispectrum 
contribution to the covariance matrix. In ^IIIl we calcu- 
late the impact of the trispectrum on the cosmological 
parameters for current and future experiments consider- 
ing a set of parameter configurations. The fiducial cos- 
mologies considered can be found in Table H] In HlVl we 
conclude with a summary. 



II. ANALYSIS 
A. Lensing Trispectrum 

We begin with a review of how to calculate the various 
contributions of the trispectrum to the covariance matrix. 
For a comprehensive introduction and derivation of these 
equations, see [14 1. 

In the flat sky approximation, the power spectrum and 
trispectrum are defined via: 



(e*(ii)e*(i2)) = (27r)2fc(ii2)cf , 

(e*(li)...e*(l4))^ = (27r)25D(ll234)T®(ll,l2,l3,l4 



(1) 



In computing these quantities the band powers are ob- 
tained, as 



subscript s^ stands for the i*^ shell in multipole space over 
which we are integrating. The signal covariance matrix 
Cij is given by 
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^'^^T«(l„-l„l„-l,)(,4) 



(27r)2 



where A is the survey area in steradians. The first term 
of the covariance matrix Cij is the Gaussian contribution 
to the sample variance and includes, in addition to the 
primary component, contributions through lensing and 
secondary effects. The second term Tij represents the 
non-Gaussian contributions contained in the trispectrum. 

The surveys we consider also include instrumental 
noise. The noise is incorporated into our Gaussian vari- 
ance as an additional contribution to the power spectrum 



(5) 



where Ni is the power spectrum of the detector and 
other sources of noise introduced by the experiment. 
These noise contributions are included in respective ex- 
periment's publicly released data and significantly im- 
pact the covariance matrix, in particular on small angular 
scales. 

For the power spectrum covariance, we are interested 
in the case where I2 — —li with |li| = Z^, and I4 = — 13 
with II3I = Ij. These conditions render parallelograms 
for the trispectrum configuration in Fourier space. 

To account for the lensing contribution to the trispec- 
trum, we must compute 



rQ(i„-i„i„-i,) = 
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~A~2n 



e(i)e(-i). 



(2) 



where As- = As{li) = J d^lWi is the area of 2D shell in 
multipole space, Wi is the i**^ window function and the 



where Cf. is the TT power spectrum and Cf is the lens- 
ing power spectrum. This equation includes all terms 
with no additional permutations. Similarly, for the 
Icnsing-secondary trispectrum effects, such as lensing- 
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ISW and lensing-SZ, wc calculate 
T«(l„-l„l„-l,) = 



2(l,.l, 



mass function of dark matter halos as rendered by [l8| 
utilizing the linear transfer function of [l9| . and y is the 
dimensionlcss two-dimensional Fourier transform of the 
projected Com pto n y-parameter, given via the Limber 
approximation [20[ by: 



(7) 



where the s is a place holder denoting either the ISW 
or SZ contribution. 



B. SZ Trispectrum 

In addition to the lensing contributions to the trispec- 
trum above, we consider contributions from the inverse 
Compton scattering of the CMB photons. The SZ con- 
tribution to the trispectrum is given by HBl : 



47rro 
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2 / xSin(;x//s) 
dxx yw[x) — — — , 
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(9) 



where the scaled radius x = r/r^ and 1^ — AaI'^s such 
that (La is the angular diameter distance and is the 
scale radius of the three-dimensional radial profile 
of the Compton y-parameter. This profile is a function 
of the gas density and temperature profiles as modeled 
in pH . Hence, we incorporate the contributions obtained 
from the SZ effect along with those from lensing, lensing- 
ISW, and lensing-SZ effects to the covariance matrix in 
Eqn. 3. 
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where gi, is the spectral function of the SZ effect, 
V{z) is the comoving volume of the universe integrated 



to a redshift of Zmax = 4, 

that [logio(Minin),logio(A/„ 



M is the virial mass such 
ax)] = [11, 16], dn/dM is the 



C. The Weak Lensing Scaling Parameter Al 



To first order in ^, the weak lensing of the CMB 
anisotropy trispectrum can be expressed as the con- 
volution of the power spectrum of the unlensed tem- 
perature Ci and that of the weak lensing potential 
C/'^'^ [H ii, [H]. The magnitude of the lensing poten- 
tial power spectrum can be parameterized by the scaling 
parameter A^, defined as 
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FIG. 1: The impact of varying the lensing scaling parameter 
on the lensed CMB temperature power spectrum, for Al ~ 
[0,2,5,10]. 



Cf* ^ AlC^^. (10) 



Thus, Al is a. measure of the degree to which the ex- 
pected amount of lensing appears in the CMB, such that 
a theory with A^ = is devoid of lensing, while A^ = 1 
renders a theory with the canonical amount of lensing. 
Any inconsistency with unity represents an unexpected 
amount of lensing that needs to be explained with new 
physics, such as dark energy or modified gravity [isl. [23|. 
The impact of this scaling parameter on the lensed CMB 
temperature power spectrum can be seen in Fig. [1] Qual- 
itatively, Al smoothes out the peaks in the power spec- 
trum and can therefore also be viewed as a smoothing 
parameter in addition to its scaling property. Given that 
Al primarily affects the temperature power spectrum on 
small angular scales, we also explore the possibility that 
it deviates from unity as secondary non-Gaussianitics are 
accounted for in the analysis. 
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ACBAR 


CBI 




72.0 


73.2 




0.02282 


0.02291 




0.1108 


0.1069 


T 


0.088 


0.086 


ris 


0.964 


0.960 




2.41 X 10"^ 


2.09 X 10"^ 



TABLE I: Fiducial model used to generate C; and C, 
CAMB for ACBAR and CBI. 
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FIG. 2: The ratio of the diagonal of the covariance matrix 
plus trispectrum effects with the diagonal of the covariance 
matrix neglecting trispectrum contributions for ACBAR (top) 
and CBI (bottom). The x in the C'p='"==+'^ is the specific con- 
tribution being added. We see that both the combined lensing 
effects and SZ alter the covariance matrix by at most a few 
percent. The choppiness is due to two issues: the covariance 
matrices provided by ACBAR and CBI are not smooth and 
the window functions used to bin the data take on different 
shapes and are themselves not smooth. 



III. RESULTS 
A. Analysis Using ACBAR and CBI Data 

We compute Cf and Cl"^, separately for each afore- 
mentioned experiment by utilizing the publicly available^ 



Fortran code CAMB [27| for the fiducial models found in 
Table m These fiducial models were taken from NASA's 
Lambda website^ which houses many best fit fiducial 
models derived from several different data sets. From 
these spectra, we use Eqns.HHHto determine the trispec- 
trum contributions from lensing, lensing-ISW, lensing- 
SZ, and SZ effects. We then add each corresponding 
trispectrum contribution to each experiment's publicly 
available covariance matrix. 

For our analysis, we use the same multipole binning 
as the aforementioned experiments, which are obtained 
from the window functions plotted in Fig. These win- 
dow functions force proper normalized bins when we in- 
tegrate over I space. To factor in the area of each survey, 
we note that the ACBAR survey covers 600 sq. degrees 
and CBI covers 143 sq. degrees. 

Fig. [2] shows the diagonal ratio of the covariance ma- 
trix together with non-Gaussian contributions, with the 
purely Gaussian covariance matrix, for the ACBAR and 
CBI surveys. We plot only the diagonal, as these en- 
tries are several orders of magnitude larger than the off- 
diagonal pieces. All contributions along the diagonal of 
the covariance leave corrections no larger than on the 
order of 1%. 

Cross-correlations become significant when the trispec- 
trum is considered, since a fully Gaussian theory should 
not have any cross-correlation terms in the covariance 
of the temperature power spectrum. As seen in Fig. |31 
the non-Gaussian corrections to the cross-correlations 
lie on the order of a factor of 100 as they correct off- 
diagonal terms that otherwise are negligible. However, as 
will be shown, despite the large corrections to the cross- 
correlations, the trispectrum effects on the cosmological 
parameters are negligible. 

Having obtained corrections to the binned covariance 
matrix from trispectrum effects, we now compute the im- 
pact of these on the cosmological parameter constraints. 
The method we use in our analysis is based on the 
publicly available Markov Chain Monte Carlo (MCMC) 
package CosmoMC [1^ with a convergence diagnostic 
done through the Gelman and Rubin statistics. We sam- 
ple the following seven-dimensional set of cosmological 
parameters, adopting flat priors on them: the baryon 
and cold dark matter densities ilf, and fic, the ratio of 
the sound horizon to the angular diameter distance at the 
decoupling, 9s, the scalar spectral index Us, the overall 
normalization of the spectrum Ag a,t k = 0.002 Mpc~^, 
the optical depth to reionization, r, and the weak lensing 
parameter A^. Wc use the same window functions and 
band powers supplied from the considered experiments 
(see Fig. U). We only use the band powers for I > 900 
for both ACBAR and CBI in our CosmoMC analysis, 
and use WMAP data for calculations on larger angular 
scales. The non-Gaussian contributions to the trispec- 



^ |http://camb.info7] 



^ http://lambda.gsfc.nasa.gov/ 
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ACBAR Normalized Trispectrum+SZ Contribution 




CBI Normalized Trispectrum+SZ Contribution 




FIG. 3: The absolute values of the Gaussian covariance plus 
trispectrum matrix divided by the Gaussian covariance matrix 
alone. In the top figure we show corrections to the ACBAR 
data set, and in the lower figure to the CBI data set. The 
trispectrum contribution to the cross-correlations is about 100 
times larger than the Gaussian contribution. 



trum are incorporated into the binned covariance matrix 
for use in CosnioMC. 

Table |TT] shows the extent to which the trispectrum 
contributions affect the constraints on cosmological pa- 
rameters for ACBAR and CBI data sets combined with 
WMAP five-year data. The trispectrum has a negligible 
impact on the cosmological parameters for these simu- 
lations. For the cases where the weak lensing scaling 
parameter and the running of the spectral index are in- 
cluded the constraints on the cosmological parameters 
change by less than 10%. 



B. The Weak Lensing Scaling Parameter 

Using a prior ACBAR dataset, it was found that the 
weak lensing scaling parameter Al is inconsistent with 
the expected value of unity at 2.5(t [15|. This led to 
an apparent revision by the ACBAR team and result- 
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FIG. 4: The window functions used for our analysis of the 
ACBAR and CBI data sets. 
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FIG. 5: The multipole-dependent ratio Ci/ACi for the con- 
sidered experiments (ACBAR in dot-dashed, CBI in dotted, 
Planck in dashed, and CMBPol in solid). 



ing in a recalibration of the data. This apparently 
removed the result first demonstrated by Calabrese et 
al. [l^ and the ACBAR team's own revised analysis 
showed Al to be consistent with unity at a significance 
around l.lcr [l6j . However, given that this parameter 
is strongly sensitive to small angular scales where other 
secondary anisotropics contribute and non-Gaussianities 
become important, we have also chosen to include this 
parameter into our analysis. We have explicitly carried 
out two runs with this scaling parameter. For the first 
run, we incorporated Aj^ without any assumption of pri- 
ors or the addition of other parameters, as shown in Ta- 
ble [III We find Al to be consistent with unity at 1.40-, 
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TABLE II: The first column represents the cosmological parameters considered. In the second column we quantify the 
parameter constraints without non-Gaussian effects. The cosmological parameters determined from the combined ACBAR, 
CBI and WMAP five-year CMB data sets along with the non-Gaussian contributions from the trispectrum are shown in the 
third column. The same data without the trispectrum but with the variation of Al is shown in the fifth column. The columns 
denoted Act represent the percent change of the la errror bars. 



Parameter 


CMBData 


CMBData±Tri. 


Aa/a 


CMBData±ylL 


Aa/a 




0.02284 ± 0.00058 


0.02283 ± 0.00059 


2% 


0.02300 ± 0.00063 


9% 




0.1103 ±0.0060 


0.1105 ±0.0061 


2% 


0.1082 ±0.0061 


2% 


e 


1.0408 ± 0.0026 


1.0410 ± 0.0028 


8% 


1.0413 ±0.0028 


8% 


T 


0.090 ±0.017 


0.089 ±0.017 


0% 


0.090 ±0.017 


0% 


Us 


0.966 ±0.014 


0.966 ±0.015 


7% 


0.968 ±0.015 


7% 


Al 


1 


1 




2.32 ±0.93 




log[10i°^.l 


3.180 ±0.047 


3.178 ±0.048 


2% 


3.166 ±0.048 


2% 




0.741 ± 0.029 


0.740 ± 0.029 


0% 


0.752 ±0.030 


3% 


Age/GYr 


13.67 ±0.13 


13.67 ±0.14 


8% 


13.63 ±0.14 


8% 




0.259 ± 0.029 


0.259 ±0.029 


0% 


0.247 ±0.029 


0% 




0.815 ±0.031 


0.801 ± 0.033 


7% 


0.789 ±0.035 


13% 




10.5 ± 1.4 


10.4 ± 1.4 


0% 


10.4 ± 1.4 


0% 


Ho 


71.9 ±2.6 


71.8 ± 2.7 


4% 


73.1 ± 2.8 


8% 



TABLE III: In the second column, we quantify the cosmological parameter constraints incorporating the additional free 
parameters Al and dus/dlnk. In the third column, we show the cosmological parameters determined from the combined 
ACBAR, CBI and WMAP 5-year CMB datasets with contributions from the trispectrum, along with the additional free 
parameters Al and dUs/dlnk. 



Parameter 


CMBData± ±Run 


CMBData±AL±Run±Tri. 


Aa/a 




0.02283 ± 0.00080 


0.02269 ± 0.00081 


2% 




0.1135 ±0.0089 


0.1136 ± 0.0090 


1% 


e 


1.0410 ±0.0029 


1.0410 ± 0.0030 


3% 


T 


0.090 ±0.017 


0.090 ±0.018 


6% 


ris 


0.997 ± 0.048 


1.001 ±0.049 


2% 


Al 


2.10 ±0.90 


2.02 ±0.94 


4% 


dUs/dlnk 


-0.014 ±0.024 


-0.017 ±0.025 


4% 


logllQ^^As] 


3.167 ±0.049 


3.167 ±0.049 


0% 




0.726 ± 0.047 


0.723 ± 0.048 


2% 


Age/GYr 


13.67 ±0.16 


13.69 ±0.17 


6% 




0.275 ±0.029 


0.277 ±0.029 


0% 


as 


0.810 ±0.039 


0.807 ±0.041 


5% 




10.8 ± 1.5 


10.8 ± 1.5 


0% 


Ho 


71.0 ±4.0 


70.7 ±4.1 


3% 



for the same reasons that the ACBAR team found this 
parameter to be consistent with the theoretical value in 
the fiducial cosmological model. 

For the second simulation, in addition to A^^ wc in- 
corporated secondary non-Gaussianities and allowed the 
spectral index to run as this parameter is sensitive to 
physics on small scales. However, the constraint on the 
lensing amplitude remains unaffected. Nevertheless, it 
is possible that non-Gaussian effects are partly respon- 
sible for the inconsistency of the scaling parameter with 
unity, and future data sets with less noise on small an- 
gular scales are required to determine this. 



C. Mock Runs 

In order to understand the extent to which the trispec- 
trum impacts the cosmological parameters in a prospec- 
tive survey with less noise, we created TT mock datasets 
assuming the best fit WMAP5 cosmological parameters, 
with the optical depth fixed at r = 0.087, and with noise 
properties consistent with Planck's 143 GHz channel. In- 
stead of considering cross-correlations, we only use the 
diagonal entries of the binned covariance matrix. 

Table IIVI shows the results of the Planck mock simu- 
lation. Since Planck has significantly less noise at high 
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TABLE IV: In the second column, we show the cosmological parameter constraints from mock Planck data. In the third 
column, we quantify the non-Gaussianity induced shift of the best fit parameter values defined as the ratio of the parameter 
shift with the uncertainty. Simulation results combining the trispectrum with the Planck mock are shown in the third column. 



Parameter 


Planck 


Planck + Tri. 


Api/a 


Act 




0.02281 ± 0.00016 


0.02271 ± 0.00018 


-0.62 


13% 




0.1122 ± 0.0016 


0.1132 ± 0.0018 


0.63 


12% 


e 


1.04095 ± 0.00037 


1.04084 ± 0.00044 


-0.30 


19% 


Us 


0.9570 ± 0.0042 


0.9538 ± 0.0047 


-0.75 


12% 


log[10iM.] 


3.211 ±0.017 


3.224 ± 0.020 


0.73 


18% 




0.7330 ± 0.0086 


0.727 ±0.010 


-0.67 


14% 


Age/GYr 


13.686 ± 0.028 


13.704 ± 0.033 


0.65 


18% 




0.2670 ± 0.0086 


0.273 ±0.010 


0.67 


14% 




0.8071 ± 0.0068 


0.8113 ±0.0076 


-0.08 


12% 


Zre 


10.395 ±0.071 


10.447 ± 0.081 


0.74 


14% 


Ho 


71.12 ±0.76 


70.62 ± 0.86 


-0.67 


13% 



TABLE V: In the second column, we show the cosmological parameter constraints from mock CMBPol data. Simulation 
results combining the trispectrum with the CMBPol mock are shown in the third column. 



Parameter 


CMBPol 


CMBPol ± Tri. 


Api/a 


Act 




0.02292 ± 0.00013 


0.02279 ± 0.00016 


-1.01 


23% 




0.1111 ±0.0014 


0.1126 ±0.0017 


1.09 


21% 


e 


1.04111 ±0.00028 


1.04110 ±0.00036 


-0.52 


29% 


ris 


0.9603 ± 0.0036 


0.9558 ± 0.0044 


-1.23 


22% 




3.198 ±0.015 


3.216 ±0.018 


1.19 


20% 




0.7392 ± 0.0072 


0.7309 ± 0.0091 


-1.15 


26% 


Age/GYr 


13.663 ± 0.023 


13.689 ± 0.030 


1.11 


30% 




0.2608 ± 0.0072 


0.2691 ±0.0091 


1.09 


26% 




0.8026 ± 0.0057 


0.8088 ± 0.0071 


-0.08 


25% 




10.336 ± 0.059 


10.410 ± 0.073 


1.25 


24% 


Ho 


71.69 ±0.65 


70.95 ± 0.79 


-1.15 


22% 



I, the trispectrum effects on the cosmological parame- 
ters is more apparent and the confidence regions of the 
most extreme changes are seen in Fig. [6l Here, the la 
parameter constraints weaken by up to 20%. Thus, the 
trispectrum has a more significant impact for measuring 
the constraints on the cosmological parameters in exper- 
iments with a noise-level comparable to that of Planck 
on small angular scales. 

Furthermore, in anticipation of CMBPol [s^l, which 
should have noise an order of magnitude less than Planck, 
we have carried out a mock CMBPol-like run with noise 
estimates consistent with the EPIC concept mission [U . 
The results of this run are seen in Table |Vl For sur- 
veys with noise level comparable to that of CMBPol, 
non-Gaussian contributions in the trispectrum weaken 
the parameter constraints by up to 30%. 

We have included in Table IIVI and Table |V] the param- 
eter shifts induced in one realization of the mock CMB 
spectrum. Although the size of the parameter shifts ex- 
tend up to 1.5(7, the shifts are not statistically significant 
as errors are correlated. We expect that for a set of real- 



izations the shifts will cancel on average. 



IV. CONCLUSIONS 

Secondary non-Gaussianities in the trispectrum from 
weak lensing, the ISW effect, and the SZ effect have 
the potential to weaken the constraints on the cosmo- 
logical parameters in experiments with high enough sen- 
sitivity to physics on small angular scales. While re- 
cent experiments such as ACBAR and CBI are too noisy 
for the trispectrum to become important at large multi- 
poles, future surveys need to incorporate the trispectrum 
to avoid over-confident parameter estimation. The un- 
derlying cosmology constrained by Planck suffers a 20% 
degradation when these secondary contributions to the 
temperature anisotropics are properly incorporated, and 
this rises to 30% in CMBPol. Accounting for secondary 
non-Gaussianities in our analysis shows the weak lensing 
scaling parameter to be consistent with unity at 1.5a. 
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